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1. The various theories of thin elastic shells which have hitherto been proposed have 
been discussed by Mr. Love"^ in a recent memoir, and it appears that most, if not all of 
them, depend upon the assumption that the three stresses which are usually denoted 
by E, S, T are zero : but, as I have recently pointed out,t a very cursory examination 
of the subject is sufficient to show that this assumption cannot be rigorously true. It 
can, however, be proved that, when the external surfaces of a plane plate are not sub- 
jected to pressure or tangential stress, these stresses depend upon quantities propor- 
tional to the square of the thickness, and whenever this is the case they may be treated 
as zero in calculating the expression for the potential energy due to strain, because 
they give rise to terms proportional to the fifth power of the thickness, which may be 
neglected, since it is usually u^nnecessary to retain powers of the thickness higher 
than the cube. It will also, in the present paper, be shown by an indirect method 
that a similar proposition is true in the case of cylindrical and spherical shells, and, 
therefore, the fundamental hypothesis upon which Mr. Love has based his theory, 
although unsatisfactory as an assumption, leads to correct results. A general expression 
for the potential energy due to strain in curvilinear coordinates has also been obtained 
by Mr. Love, and the equations of motion and the boundary conditions have been 
deduced therefrom by means of the Principle of Virtual Work, and if this expression 
and the equations to which it leads were correct, it would be unnecessary to propose 
a fresh theory of thin shells ; but although those portions of Mr. Love's results which 
depend upon the thickness of the shell are undoubtedly correct, yet, for reasons which 
will be more fully stated hereafter, I am of opinion that the terms which depend upon 
the cube of thickness are not strictly accurate, inasmuch as he has omitted to take 
into account several terms of this order, both in the expression for the potential 
energy and elsewhere. His preliminary analysis is also of an exceedingly complicated 
character. 

2. Throughout the present paper the notation of Thomson and Tatt's " Natural 
Philosophy " will be employed for stresses and elastic constants, but, for the purpose 

* *PM1. Trans./ A, 1888, p. 491. 

t ' London Matb. Soc. Proc.,' vol. 21, p. 33, 
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of facilitating comparison, Mr. Lovers notation will be employed for strains and 
directions. It will also be convenient to denote the values of the various quantities 
involved, at a point P on the middle surface of the shell by unaccented letters ; and 
the values of the same quantities at a point P' on the normal at P, whose distance 
from P is h\ by accented letters. The radius of the shell will also be denoted by a, 
and its thickness by 2h, 

The theory which it is proposed to develop for cylindrical and spherical shells 
is identical, except in matters of detail, with the theory of plane plates which I 
recently communicated to the London Mathematical Society,"^ but for the sake of 
completeness a short outline will be given. 




In the figure let OADB be a small curvilinear rectangle described on the middle 
surface of the shell, of which the sides are lines of curvature ; and let us consider 
a small element of the shell bounded by the external surfaces, and the four planes 
passing through the sides of this rectangle, which are perpendicular to the middle 
surface. 

The resultant stresses per unit of length which act upon the element, and which are 
due to the action of contiguous portions of the shell, are completely specified by the 
following quantities ; viz., across the section AD, 

T^ = a tension across AD parallel to OA, 

Mg = a tangential shearing stress along AD, 

N2 = a normal shearing stress parallel to 00, 

G3 = a fl.exural couple from to A, whose axis is parallel to AD, 

H^ = a torsional couple from B to 0, whose axis is parallel to OA. 

Similarly the resultant stresses per unit of length which act across the section BD are, 

T^ = a tension across BD parallel to OB, 
M3 = a. tangential shearing stress along BD, 
N^ = a normal shearing stress parallel to 00, 
Gj = a flexural couple from B to 0, whose axis is parallel to BD, 
Hg = a torsional couple from to A, whose axis is parallel to OB, 

* ' London Matli. Soc. Proc./ vol. 21, p. 33. 
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If the edges AD, BD were of finite length, there would also be a couple whose axis 
is parallel to the normal, but since this couple is proportional to the cube of the edge, 
it vanishes in comparison with the other stresses when the rectangle OADB is 
indefinitely diminished. 

We shall denote the components of the bodily forces per unit of mass in the 
directions OA^ OB, 00 by X, Y, Z ; but for reasons which will be more fully explained 
hereafter, we shall suppose that these forces arise solely from external causes, such as 
gravity and the like. All forces arising from pressures or tangential stresses applied 
to the surface of the shell will be expressly excluded. 

The first step is to write down the equations of motion of an element of the shell 
in terms of the sectional stresses,* which can be done by the usual methods ; we shall 
thus obtain six equations, three of which are formed by resolving the forces parallel 
to OA, OB, 00, and three more by taking moments about these lines. 

These equations will not, however, enable us to solve any statical or dynamical 
problems ; in order to do this we require the equations of motion in terms of the dis- 
placement of a point on the middle surface and their space variations with respect to 
the coordinates of that point. 

3. The values at P' of all the quantities with which we are concerned are functions 
of the position of P^, and are, therefore, functions of {r, z, (f)) or (r, 0, (}>), according 
as the shell is cylindrical or spherical. If, therefore, €^' be the value of any such 
quantity at P', and ^ the value of the same quantity at the point P, which is the 
projection of P' on the middle surface, it follows that 

= ca + .f^) + i..(^^)+ (,) 

by Taylor's theorem, where the brackets are employed, as will be done throughout 
this paper, to denote the values of the differential coefiicients at the middle surface 
where r = a. 

Objections have been raised by Saint-Venant and endorsed by Mr. Love, to the 
supposition that the first few terms of the expansion by Taylor's theorem of the 
quantities involved may be taken as a sufiicient approximation. If, however, this 
objection were valid, it would appear to me to upset the greater part of most physical 
investigations ; inasmuch as it is always assumed as a general principle, that when a 
quantity is known to be a function of the position of a point P, its value at a neigh- 
bouring point P' may be obtained by Taylor's theorem, unless some physical 
discontinuity exists in passing from P to P'. If, therefore, we put <E = R, we may 
write 

R' = A + Aj// + iA^h'^ + . . . ...... (2) 

where the A's are functions of the position of P and also of the thickness of the shell. 

* See Besant '' On the Equilibrium of a Bent Lamina," ' Quart. Journ. Math./ yoI. 4, p. 12. 
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A question which is of fundamental importance in the theory now arises, as to the 
way in which the A's depend upon h. 

If W were of the order of the square of the thickness, it is evident that A and A^ 
could not contain any powers of h lower than the second and first respectively, whilst 
Ag could not contain any negative power of h. The A's are entirely unknown 
quantities, and as there appears to be no possibility of determining them by an a 
priori method, it seems hopeless to attempt to construct any theory of thin shells 
without the aid of some assumption which will enable us to get rid of them. If, how- 
ever, we assume, as has been practically done by previous luriters^ that, ivhen the 
surfaces of the shell are not subjected to any surface forces such as pressures or 
tangential stresses, R' and also S' and T', so far as they depend on h and h', are 
capable of being expressed in thefornfi 

ivhere ic^ is a homogeneous n-tic function of h and h\ the problem can be completely 
solved without attempting to deterraine by any a priori method the values of any 
unknown quantities, and upon this fundamentcd hypothesis the theory of the present 
paper will be based, 

. There is some direct evidence of the truth of this hypothesis. In the case of a 
plane plate of infinite extent^ it can be proved to be true by means of the general 
equations of motion of an elastic solid ;^ and for the purpose of testing the hypothesis 
in the case of a curved shell, I have recently investigated to a second approximation, 
so as to obtain the term in h^, the period of the radial vibrations of an indefinitely 
long cylindrical shell, by means of the general equations, and also by means of the 
theory of thin shells, and both results agree.t But far the most conclusive evidence 
in favour of the truth of the hypothesis is furnished by the results to which it leads ; 
and I have, therefore, conducted the following investigation in such a manner as to 
furnish a test of the correctness of the final resiilts, and consequently of the funda- 
mental hypothesis by means of which they are deduced. 

Having obtained the equations of motion of a cyhndrical and a spherical shell in 
terms of the sectional stresses, all these stresses are then calculated by a direct method, 
wit-h the exception of the tensions Ti, T^, which cannot be calculated directly, 
since they involve the unknown quantities A and A^. After that the potential energy 
and the other constitutents of the variational equation are calculated, and the variation 
worked out by the usual methods. The final result, as is always the case in such 
investigations, consists of a line integral and a surface integral, the former of which 
determines the values of the sectional stresses in terms of the displacements, and the 
latter of which determines in the same form the three equations of motion. Now, if 
the work and the fundamental hypothesis upon which the theory is based are correct, 

* Lord RatleiCtH, * London Math. Soc. Proc.,' toL 20, p. 225 ; see also foL 21, p. 33. 
t^ ' London Math. Soc. Proc.,' voL 21, p. e53. 
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the variational equation will give the correct values of the tensions T^ Tq, which are 
unknown, and will also reproduce the values of the other stresses which have been 
obtained directly. This is the first test. The second test is furnished by the 
consideration that, if we substitute the values of the sectional stresses which we have 
obtained from the variational equation, in the first three of our original equations of 
motion in terms of these stresses, we ought to reproduce the equations of motion in 
terms of the displacements, which have been obtained from the variational equation. 
This is found to be the case both when the shell is cylindrical and when it is spherical; 
and I therefore think that the fundamental hypothesis is sufficiently established. 
Having obtained the values of the sectional stresses, the boundary conditions can be 
deduced by means of Stokes' theorem, which enables us to prove that it is possible 
to apply a certain distribution of stress to the edge of a thin shell, without producing 
any alteration in the potential energy due to strain. 

The fundamental hypothesis that R', S', T' may be treated as zero is not true when 
the surfaces of the shell are subjected to external pressures or tangential stresses ; for 
if the convex and concave surfaces of the shell were subjected to pressures n^, II^, the 
value of R' as we pass through the substance of the shell from its exterior to its 
interior surface, must vary from — 11^ to — U^, and consequently (excepting in very 
special cases) R will contain a term independent of the thickness. Hence the theory 
developed in the present paper is not applicable to problems relating to the collapse 
of boiler flues, or to the communication of the vibrations of a vibrating body to the 
atmosphere. In order to obtain a theory v/hich would enable such questions to be 
mathematically investigated, it would be necessary to find the values of the additional 
terms in the variational equation of motion, which depend upon the external 
pressures ; and this is a problem which awaits solution. 

It will be convenient briefly to state the notation employed. 

In the case of a cylindrical shell, OA is measured along a generating line, and OB 
along a circular section. In the case of a spherical shell, OA is measured along a 
meridian, and OB along a parallel of latitude. 

The three extensional strains along OA, OB, 00 are denoted by cti, cr^, org ; and the 
three shearing strains about those lines by ijt^, tat^, xn-g. We shall also use the letters 
X, jm, j9 ; X\ [x, p' to denote the first and second diflerential coefficients of or^, org, ct^ 
with respect to r, when r = a. We shall also write 

E = (m — n)/{m + ^?), . K = o-i^ + o^s? 

= OTi + E (o-i + O-g), 3$ = OTg + E (o-i + (7a), 

= X + E (X + ^), ^ = ;. + E (X + ^), 

3E'= V + E (V + /), ;(?' zz: ^' + E (X' + fi). 
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Cylindrical Shells, 

4. Before we can obtain the equations of motion or tlie potential energy, it will be 
necessary to ascertain the values of the first and second differential coefficients of the 
displacements with respect to r when r = a. We shall, therefore, proceed to calculate 
these quantities. 

Putting X, /x ; V, //.' foi* the values of {dcri/4r)^ (dcr^/di^) ; (c^^o-j/cir^), (d^cr^/dr^) 
when r = a, we have 

B' = (m 4- n) cr'g + (m — n) {cr\ + (r\^ 
== (m + ^^) 0^3 + (^^ — ^^) (<^i + o-z) 

{m + ?i) f — - ] + (^ — tZ') (X + /^) f ^' 






'^2 



U/3 



But from (2), 
whence 



+ i j('^^ + ^^0 (^j + (m - n) (X' + /xO \li' + 

R^ == A + AJi + iAJi' + . . . , . 
A = (m + n) (T^ + {m — n) (cr^ + o-g) ^ 

A, = (m + ^^0 ( T^J "^ (^ "^ ^'0 (X + m) l 

A^ = (m 4" ^^) 7::i + (^'^'^ — ^) (X' 4- /^O 



* « « 



rf 



^ 



(3). 
(4), 



(5), 



where A, A^ do not contain any lower powers of h, than h^ and h respectively, and 
Ag does not contain any negative powers ofh 

If u, v\ %D be the component displacements of any point of the substance of the 
shell in the direction z^ <p, r, the equations connecting the displacements and strains 



are 



whence if 



du^ 



0-1 — 



cr 



di 



(T Q 



137 1 ' — 



1 /di/ 
T \d(p 

dio' 
dr 

di/ 
dr 

dw' 



+ v/ 



L 



v^ 1 dw' 
r r d<p 



m 



w 



dto^ 
dz dr 

1 diif dv^ 
r d(j) dz 



(6), 



taaf^ 



"P 



in — % 
m + n 



K =^ cfi^ 4* 0*^ 



6 » 



(n> 
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we obtain 



di\ 



dw 






dz 



dv 
dr 



ZZZ tSTj^ *|- 



V 



a 



1 dw 
a dcf> 



'dw" 
dr 






m -{■■ n 



- EK (8), 



and 



\ 



dr^ 

'dh' 
dr"^ 

'dhv 
'dj^ 



'd 



Wi. 



dr 

dxjs-t 
dr 



dz dr 






w. 



a 



/dm A 

\ dr j m -\- n dz 

1 ^A E c^K 



dA. ^ dK^ 



tvZ 



a{m + n) d(f> 

E + )u.) 



'a d<f> 



r 



(9). 



5. We ca,n now obtain tbe equations of motion in terms of the sectional sti'^sses. 

Let dB be an element of the middle surface whose coordinates are (a, z, ^), and 
dS' an element of a layer of the shell whose coordinates are {a + h'^ z, «^) ; then 
dS^ = (1 -}- h/jci) dS, If we considei? a small element of volume bounded by the two 
external surfaces of the shell, and the four planes passing through the sides of dS, 
which are perpendicular to the middle surface, we obtain by resolving parallel 
to OA, 

I (V 8(^) Sz + ^ (Ml 8^) S<}> t= p dB f ^ (i? - X) (1 + hr/a) dJi (10). 

But 



dip 



^h 



u 



/iu\ 



» + *' i^j + 4 ''" (£ 



accordingly if we substitute the values of {du/dr) and (d^ujdr'^) from (8) and (9) and 
recollect that all terms which vanish with h may be omitted when multiplied by h^, 
the right hand side of (10) becomes 



pc^S 2/^(i;--X) + i^'E 



^K 2%^ dw 



dz 



3a dz 



> • 



Resolving parallel to OB, 00, and then taking moments about OA, OB, 00, we 
shall obtain in a similar way five other equations, which, together with (10), may be 
written 



d^ 

dz 



J JL iA/ XVJL "I 



a d(f> 



2h {u -- X) + i h^E 



dK 2}}? dw] 



dz 



*da dz 



a d^ a dz 



p \2h(v 






6a' 



dm' 
d^^ 



d¥^ 1 d¥[^ 



dz 

1 dxJTY 



^.i-l'^"! 



2A». 



2A(w .^ Z) — |A3E(X+/x} — — EK I y (11), 



3a 






d(j:<^ 1 (i^Hg 

dz a d<}> 



3d \d(l) 

• • • • 

■y-r Q 7 Q / ail) U -A. \ 

-^,= -spf^- [^-a + aj 



^ 
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These equations will not enable us to solve the problem in hand ; in order to do 
this we require the equations of motion in terms of the displacements, and also the 
values of the sectional stresses in terms of the same quantities. 

6. The values of the couples, and also the values of M^ M^, can be obtained by 
direct calculation ; but the values of T^ T^ cannot be so obtained, since they involve 
the quantities Kh and AJi^, which are unknown, and which cannot be neglected. We 
shall, therefore, be compelled to find the expression for the potential energy, and 
employ the Calculus of Variations. 

The following results will, however, be necessary hereafter. If P', Q', R', S', T', 
U^ be the stresses at the point a + h\ z, ^, we have 



T,a S<f> = r P' {a + Ii) S<j> cW 



'k 



W 



HwlJ 



■k 



■p J^h' (^ + 1 h''^ ^"^'^ 



li! 



dr 



ce 



Ti=: 2hV + ^W 



f(PY 






+ 1 



dr^ 



2W fdV 



>{l + ~]ah(f> dh\ 



3 a \dr 






d")"^ 






m, 



di 



>i'^>4 



> 



G,= 



dr — 



til — 






2_ J^S pQ 
3 \^^ 



i¥ 



dr 



p 

^ a 









TJJ 



3 



a 



t » \ L ^ i» 



From the third, fourth, and last of these we see that (Mg — M^) a = H^, as ought 
to be the case. 
Let 

= 0-1 + E ((Ti + O-g), 3$ = CTg + E (or^ 4- cr^) 

= X + E(X + ia), ^ = ^ + E(X + ^) 

Then, in the terms multiplied by A'^ we may put 



(13). 



fdV 

\dr 






Q = 2nM, 
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whence, if _p = {dzsjdr), the last four of (12) become 



Gi = - I n¥;^, 



G3 = I nh^ 



+ 



a 



Hi=-f«/i«(i> + ^^)' TL,^lnh'p 



Y 



(14). 



Since the couples are proportional to the cube of the thickness, it follows from the 
fourth and fifth of (11), that the normal shearing stresses N^^, N3 are also proportional 
to the cube of the thickness, and therefore the terms of lowest order in the expres- 
sions for the shearing strains 77t{, 7s<l are quadratic functions of h and lfi\ since such 
functions when integrated through a section of the shell, give rise to quantities 
proportional to the cube of the thickness. This is consistent with the fundamental 
hypothesis. 

The next thing is to calculate the values of the quantities X, />c, p. 

From the first and fifth of (6) we obtain 



X = 



'di 



dr 



dxfjc^ 

dz 



dho 
dz^ 



and, since the terms in X are all multiplied by h^^ we may put 

d^w 






dz^ 



(15). 



Similarly from the second and fourth of (6) we obtain 






Lastly, 



3(S^ + ^)'~!^^i + ^^> 



(16). 



P = 



'd 



t»Tc 



d^Vj 



dr 



1 + ( - 



dv 



a \dT d<pj \dr dz 



1 du 

a^ d^ 



or 






2 d^w 1 dv 1 du 
a dzd^ a dz o? d^ 



(17) 



We have, therefore, completely determined the values of the couples in terms of 
known quantities. 

We shall also require the values oi {d^cTijdr^), {(Pcrjdr)^ (d^mjdr)^ the first two of 
which we have denoted by X', /x' ; and the last of which we shall denote by p\ The 
values of these quantities can, by a similar process, be shown to be 

d^K 



V=:E 



/x = — 



dz^ 

2fi . E d^K E ,, , , 
4" 3 T7^. — - (^ + /^) 



V 



a 



a^ d<j>^ 



a 



P 



MDCCCXC. — A., 



p V5^ 2'E d^K 

— —I— — U- — 

a Or a dzd(f) 

3 L 



(18). 
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Equation (17) and the last of (18), combined with the third and fourth of (12), 
determine the values of M^, M^. 

It will be desirable to point out at this stage of the inyestigation, that we have 
obtained materials for the complete solution of any problem in which T^^ and u are 
zero, and none of the quantities are functions of z. The boundary conditions at a 
free edge will be discussed in § 11, and the reader who does not wish to be troubled 
with the long analytical process of finding the potential energy and working out the 
variational equation of motion, may pass at once to § 10, and the following sections 
where certain problems of a fairly simple kind are discussed. 

7. We must now find the potential energy due to strain. 

By the ordinary formula, the potential energy of a portion of the shell is 

'2 



W = ifff^J(m + ^)A'^ 



+ n {m^^ + t^/^ + z^P - 4 (cr/a-/ + (x^cr^ + 0-3 V)}] (1 + h'/a) dh'dS (19), 

where the integration with respect to z and (f) extends over the middle surface of 
the portion considered. In evaluating this expression we may at .once omit mi\ mc[, 
for since they are quadratic functions of h and h\ they will give rise to terms which 
are proportional to h^^ which are to be neglected. 
Since 

it follows that 



A'=A+h'^~^^+ih'^c^,)+... 



from which it is seen that W is expressible in a series of odd powders of h. 
From (5) we obtain 

A = (1 -^ E) (oTi + <T,) + ^ 



m + n 



1'^] = (1 _ E) (\ + a) + -^ 

\ drj ^ ^ ^ ' ^ m + n 

'#A\ ,. _x ,. , . .X . A 



= (1 - E) (V + i^O + 



3 



dr^ J ^ ^ ^ ' m -\- n 

and, therefore, the portion of W per unit of area of the middle surface, which depends 
upon A', is 
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4# 



m -i- % 



h (o-i + o-^ + A/2n)^ + ^h^\+ IJ. + Aj2nf 



2¥ 



+ i^^ (^1 + 0-3 + A/2n) (X' + / + A^/2n) + -^{c^ + (t^ + A/2n) (X + /^ + Aj2n) M 20) 



So- 



in which in the last three terms we may omit the A's since they are multiplied by h^. 
Again 



whence 



2n 



'h 



-A 



47iA^ 



+ "^T" (^^3 + /^^i) (21) 



3a- 



Also 



(o-/ + or/) 0-3' (1 + h'/a) = (0-1 + cr^) erg + /^'^ (\ + ;.) (^) + l/j'^ (X' + /^') (Tg 

(X + //,) org + (cTi + (Tg) f -^ j 



+ i ^'M-. + -.) f^) + f < 



whence 



2n\ (cr/ + ^2) ^3' (1 + h! ja) dK = 4ni^ {cti + cTg) 



A 



m -{- n 



E (<Tj -f- (Tc^ 



Snh^ 



-tnhm{X + iif-inhm{a'^ + a,){\' + ix')^-~-E{X + iJL){o'^ + a,) (22). 



Lastly 



da 



|n cr3'^ (1 4- h'/a) dh' = n^ts-g^ + -| n¥p^ + f nh^m^p' + —^ts-gp 



-A 



3«^ 



(23). 



Substituting from (20), (21), (22), (23) in (19), it will be found that the term Kh, 
which is (or at any rate may be) proportional to h^^ disappears ; and thus the value 
of the potential energy per unit of the area of the middle surface is 



+ I nW {X^ + ^3 ^ E (X + i^)' + i p'} 

a 



+ f^^ (m+23/x+ltnrp) 



(24), 



in which m is written for tzrg, the suffix being no longer required. 

This is the expression for the potential energy as far as the term involving h^. The 
first term depends solely upon the extension of the middle surface ; the second term 
depends principally upon the quantities by which the bending is specified, and the 
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third and fourth consist of the products of the extensions and the quantities which 
principally depend upon the bending. 

8. Having obtained the value of the potential energy, we must in the next place 
form the variational equation of motion. This equation may symbolically be written 

SW + S€; = SU + SE . . . . . . . . (25), 

where 8W is the term which depends upon the time variations of the displacements, 
SU is the work done by the bodily forces, and S% represents the work done upon the 
edges of the portion of the shell considered, in producing the displacements, St^, Si;, Sw^ 
by the forces arising from the action of contiguous portions of the shell. It, therefore, 
follows that Sit is a line integral taken round the edge of the portion of the shell 
which is being considered ; and as one of our objects is to calculate the values of the 
sectional stresses in terms of the displacements by means of (25), it will be convenient 
to apply the variational equation to a curvilinear rectangle bounded by four lines of 
curvature. 

We must now calculate SC. We have 



'k 



Now 



SE = /) J j (u' Su' + v' Sv' + w Sio') (1 + h'/a) dN dB. 



u' 8ii (1 + h'/a) dN = 2hu hu + %h^ ~- ~^— + 3" ^^ ( ^ TT" + ^-^ S11 



dr dr 



dr^ dr^ 



2h^ /•• dBu du ^ 

I. I.. ■ I oj - » Iiiw O?/ 

Za \ dr dr 



^ 7 cs . o 7 Q dvj dSw , , 7 Q-T-^ /* • dBK. , dK. cn 



dz dk 



dz 



dz 



2K^ /• • dhiu , dw ^ 

•— — j Vj -4- ■ — ou 

Za \ dz dz 

by (8) and (9). Treating the other terms in a similar way, we shall find that the 
value of h% is 



SC = 2pA {to hu-^- vhv -\- ID Siv) dB 



+ iph'\ 



f* r 



• • • • 

dw dSw 1 fdw 
dz dz a^ \d(fi 



— V 



'dSw 
d4> 



8y +E^K8KMS 



+ \ phm J J j 1^ -^ + - -— - ^2^ (SX + S/^) + ^ 8m -f - — St^ 



dz 



a d^ 



» • . • 
(X + /x) Zw 



d^ 



2pW 
3a . 



'*'"-- d^w . dw ^ , V fdSw CN \ . 1- /^'^ \ CO 



u 



dz dz 



_ / __ — ^-„-., ^/jj I — j— — / — 
a\d(j) J a \d^ 



+ E(ti?SK + kS?^)lc^S • . (26). 
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We must, in the next place, calculate 8li. We have 



8a= {P' Su' + W 8v') {d + h') dh' d<)> + (Q' Sv' + V Su') dh' dz 

+ \l^^Bw/ad^+ [l^iBtvdz • (27). 

From the way in which SC has been calculated, we see from (8), (9), and (12), that 



f P' Su' (1 + h'/a) dh/ =T^Su+G,^^ + ^ ¥F ^'^"^ 



dr 

dSw 



dr^ 



= T,Bu-G,"^ + inhmm 



dz 



dSK 
dz 



Treating the other terms in a similar way, we find 






T,St. + M,8t.+ N,Si.-G,^+^('^ 



dz 



w 






dz 



4" H Mj St^ + T^ St; + N^ Stc? 4" 



Gj fdSw ^ 
a \d^ 

2n¥ 



TT 



% 



da 

dSw 
dz 



d(^ 



ad^ 






Za 



E3$ -.TT-' + i nW&nr 



a^ 



dz 



>dz 



• * 



(28). 



Lastly, since the shell is supposed to be so thin that X, Y, Z, may be treated as 
constants during the integration with respect to h' , 



8U 



j r (X W + Y 8v' + Z hw') (1 + h'la) dh' dB 



= 2oh 



(X 8m + Y Si; + Z St(;) cZS 



+ ip ¥E f f 

3a JJ 






if a: 



^ c?^ 



[■dS 



^ dBw Y /<^Sw? 
dz a \d^ 



Sv] + ZESK\dS 



. (29). 



9. We have now obtained all the materials for the complete solution of the problem, 
and we shall proceed to work out the variation in the ordinary way. 

Let us denote the four terms of the expression for W given in (24), when 
integrated over a curvilinear rectangle bounded by four lines of curvature, by W^, W^, 
W3, W4. Then 

r r 

SWj = Anh {^BcTi + MBcr^ + | tarStji-) a dz d^. 
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Substituting the values of a-^, o-g, tot from the first, second, and sixth of (6), and 
integrating by parts we shall obtain 



8Wi = ^nh [{^hu + 1 n^hv) ad^ + Inh J(-|- r;^hu + l^hv) 



4c7ih 



} l\ dz 2a d(pl \a a<p 



+ i 



2 



'^^\hv- 



d 



dz 



dz I 



Sw 



a 



adzd(j) , . (30). 



Now SWg, SWg, SW^ depend upon h^ ; if therefore we substitute in (25) the value 
of SWi from (30), and the portions of SC, SU, and S|i, which depend upon h, we shall 
obtain the approximate equations 



Ti = ink % % = inh 2$ 1 



Ml = M^ = 2nAT*T 



and 



(31), 



pu 



2nl'^ 



\dz '^ 2a d<f>j "^ P 



t& d(f> ^ dzj 



pv = 

pw =z ^ 2n ^ja + pZ 



+ pY 



^ 



(32). 



These equations are the same as those obtained by Mr. Love,"^ and which are 

employed by him in. discussing the vibrations of a cylindrical shell. The complete 

• • . • • • 

equations giving pu, pv, pw in terms of the displacements and their space variations 
contain certain additional terms involving W (since the common factor h disappears) 
which it is our object to determine ; but, since we do not retain terms higher than W, 
we may, if convenient, substitute the above approximate values in all terms of (25) 
which are multiplied by h?. 
Again 

SW^ = in¥ f [(IE SX + dF V + ^pSp)adzd<(>. 
Substituting the values of X, /x, p from (15), (16), and (17), we obtain 



IflE hXdzdcjy = ^ fflE ^ dz d4 



= \{fj^ "- ^^)^^ "" ff^8^^^# 



(33), 



also 



* ' PHI. Trans.,' A., 1888, pp. 538 and 540. Equations (32) correspond to Loye's equations (86), 
(87), and (88) ; and (31) to (101). 
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l£^ni.dzd<f>=- ^Jfd? {^^+ Sw + ^(a 



dBu d Sv 
+ 



dz 



d^ 



■^ Sw]> dzdcj) 



a 






a 



3 



I 



dz 



^«f^»+^^^^-i^+^+E^i8'" 



dzd^ (34), 



9,Ilu. 



1 

2 



pSpdzd(l> 






2a3 



1 

2a 



dzd,(f> 

1 ({(dp 



fpS^#-^fi'S^c^.+^ff(|8« 



a J- St; ) tfo c?<^ 



1 



p -r-77 rf^ rf<p. 
■^ dzd<p ^ 



The last integral can be evaluated in two different ways, according as we integrate, 
first, with respect to ^, and secondly, with respect to 2; ; or first with respect to Zy 
and secondly with respect to ^. The proper way to deal with such a term is, to 
evaluate the integral in both ways, and then multiply the two values by y8 and 1 — /3 
and add, where ^ is a quantity which must be determined from the conditions of the 
problem in hand. We shall thus find that the value of /3 is ^ ; we therefore obtain 



1 

2 



\p Sp dz d<f) 



2a 



^ . dp ^ dBia\ , 



1 f / ^ dp ^ , dBw\ , 



^-Lm^u-a'^S 



2a^ J J \d^ 



dz 



V 



2a 



d^p 
dzd(p "^^ ) 



hw\dzd<l> 



• • ♦ \.oo\* 



Collecting all the terms together from (33), (34), and (35), we finally obtain 



SW 



2 



^nhP 



I 



a \dz 2a dcp^ 



dSw 
dz 






+ f nW 



p 

2a 



Su 



a 



a \a d^ '^ dzj 



a \d^ 



Sv 



-, dSw 1 , 
f p -7- y dz 
'^■^ dz 



fgifjtmmr 



4nh? 
3a 



n 



az La acp 



\a dcp '^ dz] 






ydw 



a dz dcj) (36). 
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Let 



then 



'=X' +E{K'+iJL% ;^' = fi' + E{X' + fi') .... (37), 



SW3 = f n/i^f j"(iE'Sa-i + #'80-3 + 1 p'St^ + ftSV + MSfi' + I tjtS/) c^S ; 



from this result, together with (24), it is seen that SW3 and SW4 each consist 
of two parts, which may be denoted by SW^', SW3'' and SW/, SW/' respectively. The 
values of SW3' and SW/ may at once be written down from (30), by changing ^, 3$, 
OT into W, ^\ p' and IE, dF'? V respectively, and by altering the coefficients from 4:nh 
into %nh? and Anh^/Ba respectively. With regard to 8W3'' we have 



SW^' = I n¥ 




8X' +MSix' + ^ m8p') c/S. 



Substituting the value of X' from (18) and integrating once by parts, we obtain 






dz^ 



f d?K 



(M dSK 
dz dz 



rfb. 



Treating the other terms in a similar way, we shall finally obtain 



8W3'' 



^nhmn^^ + '^-Rj^^ 



dz 



$a 



d(j> 



) «# + {(^ E3$ ^ + i nhm.. ^]dz 



3a 



-%n¥\{\E(^:+~~Y-^-\--(-^ 



3 



d^ 



1 dm\ dSK 



Eas 



dz 



2a d(f>J dz a\ad(j) ^ dz / d(p 



a 



(8X + 8j^)Us 



+ 2. 
3 



nh' 



a 



jlj-. 23$8/.L + J (8^- a8|))|dS . 



(38). 



If in the first surface integral in this equation, we substitute the approximate values 
of the coefficients of dBK/dz^ &c., from (32), which we may do, since this integral is 
multiplied by h^, and then substitute the values of SWo'', 8C 8U, and 8% in (25), it 
will be found that all the terms involving dSK/dz, d8K/d(f), and 8X + 8/x cut out ; we 
are, therefore, no longer concerned with them, and the value of 8W3'' reduces to the 
last line. On this understanding we may, therefore, write 



SW'' + 8W. 



// 



4:nK' 



3a 



SX + J™ 8w + i' ^8p) dS 



4:a 



4:n¥C 



Ba i I dz 2a \d^ /J 



+ 



4:nM 

'~3a 



li 



dM dhw 1 dw fdhw 
dz dz 2a dz\d(^ . 



8i;)Us . . , (39). 



CYLINDRICAL AND SPHERICAL THIN ELASTIC SHELLS. 



449 



We are now in a position to test the correctness of some of our work, for picking 
out the terms involving dZwjd^ — hv, dSw/dz in the line integrals in (36) and (39), 
and equating them to the corresponding terms in the value of 8% which is given by 
(28), we see that we have reproduced the values of the couples, which we have 
already obtained in equation (14). We may therefore leave the couple terms out 
hencfeforth. 

Collecting all our results from (26), (28), (29), and (39) the variational equation 
becomes 



SWi + 8W, + SW3' +SW4; + I ^ [I 

r r 



d^ dhiD 1 dm fdhw ^ \1 ^q 
dz dz 2a dz \ d(p J J 



+ t''*11{(f-:)*+if*-2^)(*-S»)+E(EK-«,/«)SK 



a^ \d^ 






\ d(f> 



dS 



^iphm\\{§Bu + lf^Bv-{X + ,)8^.}d8 



2p¥ 



f{ 






dz 



a \d^ 



dS 



dSio , Y /dSw ^ . . --,_, ^^^ 



a \d<p 



+ j (T^Su + M^Sv + N3 §?/;) a c^</> + f (M^ 8u + 138^; + ^^Sw) dz . . (40) 



where the values of SW^^, SWg are given by (30) and (36), and the values of SWg^ 
SW/ are obtained from (30) by changing certain letters as we have explained above. 

We have now got rid of all the terms involving the second differential coefficients 
of Su^ Sv, Sw, and all that now remains to be done is to integrate by parts the terms 
which involve the first differential coefficients. Putting: 



2n d^ , dw u , X 
pa dz dz a a ' 



^ n dm , 1 Idw -A . Y 

/8 = --:,- + - bi - 2^ +" ' 

pa dz a \d^ j a 



we have 

MDCCCCX. — A. 



y == E (EK - wja + Z/a) • (41). 



3 M 
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/ /• / J 



f pA^ 



06 



d Sio B d Siv 



dz 






a d<p 



+ ySK]d8 



f ph^ j {y Su -{- a Bw) ad(j> ~\- ^ ph? (y Sv + y8 hw) dz 

'da 



■~ -I p/^« 



-r ^u + " "7 8'^' + » , 

a;;; a dcp \ d. 



a dcp a] J 



. (42). 



Substituting the values of SW^, SW^, 8W3', SW^^', and the right hand side of (42) 
in (40), and picking out the line integral terms, we obtain the following equations 
for the sectional stresses, viz., 



Ti 



M 



3 



N 



3 



4riAE - ^ EdiF + I- n¥^ + ^^'^ 



3C6 



3a 



A ^;,B ('5 + 1 '^P^ 



+ ^ E (aEK - u- + Z) 



C?^ 



2a (i^^ 



Go 



Hj := 



|- 7^A^ 



+ 



a 



f "A3 (p + ^) 



3 a 



, 47^^^ d,% , ^ ,^ /c^i6' 



3 a dz 



dz 



u X' 
a a 



•\ 



> 



(43) 



-/ 



which give the values of the sectional stresses across a circular section ; and 



M. :^ 2nh7^ + \ n¥'p 



•> 



T„ = 



Ni 



4.nhM - ^' E4? + -I n/i'd?' + ~ E (aEK -w + Z) 



3a 



1^1^' 



H 



2 



3 a 
^ \« d^ ^ 2 dzj ^ :ia dz^ 'da \d(j) 



2i; + Y 



^ . . (44) 



_x 



which give the values of the sectional stresses across a meridian. 

If we compare these equations with the third and fourth of (12), with (14), and with 
the fourth and fifth of (11), we see that we have reproduced (i.) the values of M^, M^ 
given by (12) ; (ii.) the values of the couples given by (14) ; (iii.) the values of the normal 
shearing stresses which are obtained from the fourth and fifth of (11), by substituting 
the values of the couples from (14). We have thus subjected our fundamental 
hypothesis to a fairly searching test. It is, however, in our power to subject it to a 
still further test ; for if we equate the coefiieients of ^u, Zv, 8w in the surface integrals 
in (40) and (42), we shall obtain the equations of motion in terms of the displacements. 
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and on substituting the values of the sectional stresses from (43) and (44), in the first 
three of (11), we ought to reproduce the equations of motion in terms of the displace- 
ments which we have obtained from the variational equation. 

From (30), (36), (40), and (42) it follows that these equations are 



Pi 



dK 



2 (t^ - X) + ihm p - 



dz 



2}i? div 

3a dz 



, :d^ , 1 dm 

= 47i I — + 

dz 2a d(j>^ 



Anlh? ^ d^ 

oCt CLZ 



+ f "i^h^ 



'd 



dz 



^ 2ad(pJ ^ ?,a dz ^ ^ P dz 



(45), 



p 2 (r 



Y) + ;^ E -; + 



3(X 



d^ 



dvf 
3a \ d(f)^ 



, ,1(^23 , T dm\ , 4nh^ j , 



3a2 



+ inh^(lM^^^P' 



d'^ 



^) ^1 + ^ .f + 2 ^^ 



dp 

dz 



d 



m 



\a d^ 



dz 



_L M'Ay > M' /4^' _ 9^ -1- Y 

^ 3a d<t> ^ 3a^ \(^^ "^ ^ 



(46), 



Pi 



2('mj — Z) — lh^E(i-\-'il) 



3a 



EK 






a 



I J ill « 



Sa^ 



a' 



dz' 



I 7 j9 1 l^d)' "{" <^ 



<f^J[? 



^^3 ' -^ ^ ^'^ dzd(f) 



- ^V' + *f 7(f + r3 + tM^' ft - ^ + ~ 

oa dz \ dz 2a d^/ ' ^ ^ dz \dz a a 



oa 



2ph? d fdvj •• ^ 

"^ ^^l4\d4 ~2i; + Y 



If E KK - ;« + z) . 



(47), 



If we compare these equations with the equations obtained by substitutmg 
the values of the sectional stresses in the first line of (11), it will be found that 
they agree in every respect. 

10. It will hereafter be necessary to consider certain problems inwhich the middle 
surface is supposed to experience no extension or contraction throughout the motion ; 
and it will, therefore, be necessary to obtain the requisite equations when this is 
supposed to be the case. 

The conditions of inextensibility are 



<T, 



0, cr^r=z 0^ ts- = ; 



or 



dib 
dz 



= 0, 



dv 
d(p 



3 M 2 



dii- , dv 
d(f> dz 



• (48), 
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wliich require that 



u 



Ua, 



V 



% 



^ 
d^ 



! 

"r 



d}\] dN 



w 



z 



d(f)^ d(j) 




where U and V are functions of <^ alone. 

In this case the potential energy reduces to the second line alone, and from (15), 
(16). and (17) we obtain 



X— 



/x 



"A 



1 fd^w , 
I — _[_ 10 



a" \d(l} 



r3 



3 



* > 



P 



2 / d^w 
a \dzd(p 



dv\ 

dk 



and from (24) 



W 



^nh^ 



m 



'd?"W 



Sa^ [a^ (m -\- n) \rZ^ 



3 



r^ y;\ ^ 



dhu 



dv^^ 



dzd(f> dz^ 



(50) 



• • (51), 



which agrees with the expression obtained by Lord Rayleigh.* 
Also from (14) 



Gi 



Go 



Hi 



4«-^^^(i + E)f^:+^^'' 



3a2 



^cl(f> 



3 



2 — — -| nh^p 



V 



. (52). 



-^ 



The values of the stresses M^, Mg may be obtained either from (43) and (44), or 
from (12) combined with (15), (16), (17), and (18) by introducing the conditions of 
inextensibility ; and the values of Tj, Tg might be calculated by taking the variation 
subject to the conditions of inextensibility, and using indeterminate multipliers. 
This process would not, however, be of much assistance, inasmuch as it would intrd^ 
duce two undetermined quantities into the values of T^, T^, which depend upon the 
boundary conditions; whereas in this case the values of T^, Tg can be obtained directly 
from the first and third of (1 1) combined w^ith (49). The values of N^, N^ are given 
by the fourth and fifth of (11) combined with (50) and (52). 

11. We must lastly consider the boundary conditions. 

Equations (43) and (44) determine the stresses on the line elements arf^ and ;c?^ 
respectively, which are produced by the action of contiguous portions of the shell ; 
and it might at first sight appear, as was supposed by PoissoN,t that when a shell 



* ' Roy. Soc. Proc.,' vol. 45, p. 116. 

t ^ Paris, Acad, des Sciences, Memoires,' 1829, vol. 8, p. 357. 
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of finite dimensions is under the influence of forces and couples applied to its edges, 
these equations would give the values of such forces or couples, and that the 
conditions to be satisfied at a free edge would require that each of the above five 
stresses should vanish at a free edge. Kirchhoff'''" has, however, shown that this is 
not the case, but that the boundary conditions are only four in number ; and the 
reason of this is, that it is possible to apply a certain distribution of stress to the edge 
of a shell, without producing any alteration in the potential energy. 
By Stokes' theorem, 



7iH 






, dSw\ 
d(p 



dij) + 



f/dW 



dz 



8w + H' ^~r~ ] dz = ; 

dz 1 



the integration extending roimd any curvilinear rectangle bounded by four lines of 
curvature OA, AD, DB, BA. If, therefore, we apply to the side AD the stresses 



M,' = H'M N,' = \ ^ 



Ml = H ; 



to the side DB the stresses 



1 '"" ~d^' 



2 



H', 



and to the sides BO, OA, corresponding and opposite stresses respectively, 
preceding integral becomes 



-< 



Mo' ^v -+ N/ Ziv + 



H/ Idhw 



s 



a \d(j) 



ad^ + ( N^' Sio — H^' — — ) (^2; ™ 0, 



dz j 



which shows that the work done by these stresses is zero. Such a system of stresses 
may, therefore, be a-pplied or removed without interfering with the equilibrium or 
motion of the shell 

Let us now suppose that the rectangle OADB, instead of being under the action of 
the remainder of the shell, is isolated, and that its state of strain is maintained by 
means of constraining stresses applied to its edges ; then it follows that if, instead of 
the torsional couples H|, fl^, due to the action of contiguous portions of the shell, we 
apply torsional couples l^^, ^3, where 



I, = H, + H' 



#««cee«4* 



TT 

:2 — -^3 




^ ^ Orelle,' vol. 40, p. 51, 1850, and Collected Works, p. 237. 
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the state of strain will remain unchanged, provided we apply in addition the stresses 






(55) 



and 

— — X 1 1 "~| ' -J « * . , , , , , . \tjOfj 



d:--' 



whence, eliminating H' between {5B), (55), and (54), (56) respectively, we obtain 

#l,a - 111 = M,a » H, ^ 

^ a d(j> ^ a d(f) 

and 



'^'^ dz '^'^ dz • = •'.-.. {06). 



In these equations the Roman letters denote the stresses due to the action of 
contiguous portions of the shell, whose values are given by (43) and (44), whilst 
the Old English letters denote the values of the actual stresses applied to the 
boundary. If, therefore, the shell consists of a portion of a cylinder which is bounded 
by four lines of curvature and whose edges are free, the boundary conditions along 
the circular edges are obtained by equating the right hand sides of the first and 
fourth of (43), and the right hand sides of (57) to zero, the first two of which express 
the condition that the tension perpendicular to, and the fl.exural couple about, a line 
element of the circular edge must vanish when the edge is free ; and the boundary 
conditions along the straight edge are similarly obtained by equating the right hand 
sides of the first, second, and fourth of (44), and the right hand side of (58) to zero, 
the first three of which express the conditions that the tangential shearing stress, the 
tension and the flexural couple must vanish when the free edge is a generating line. 
We may also, if we do not wish to introduce the time and the bodily forces into these 

• • • • • • 

equations, substitute for u — X, v — Y, lo — Z their approximate values from (32). 

12. We have now obtained all the materials we require, for a perfectly accurate 
approximate solution of any problem relating to the vibrations of a thin cylindrical 
shell as far as the terms involving the cube of the thickness, but before proceeding 
to discuss any problems, it will be necessary to make some remarks respecting 
Mr. Love's paper. The first line of my expression for the potential energy which is 
given in (24), and which involves h and not h% agrees with the expression obtained by 
Mr. Love and other writers ; also the approximate equations of motion (32) agree, as 
has been already pointed out, with the corresponding equations obtained by him, and 
by means of which he has discussed the extensional vibrations of a cylinder. It will 
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also hereafter appear, that observations of a precisely similar character apply to the 
corresponding equations which determine to a first approximation the extensional 
vibrations of a spherical shell. This portion of his paper therefore appears to be per- 
f(3ct]y satisfactory ; but that portion which involves the terms depending upon the 
cube of the thickness is open to criticism. 

In the first place, he appears to have employed a system of rectangular axes, con- 
sisting of the normal at a point on the middle surface, and the tangents to the two 
lines of curvature through that point. Now, although it is immaterial, so long as 
we confine our attention to infinitesimals of the first order, whether a quantity is 
measured along the tangents to three orthogonal curves or along the curves them- 
selves, yet when it is necessary to take into consideration infinitesimals of higher 
orders, which is always the case whenever an investigation involves changes of 
curvature, a method in which everything is referred to rectangular axes requires care ; 
and on comparing the terms in h^ in (24) with the corresponding terms in Mr. Love's 
expression for the potential energy, it will be seen that he has omitted several terms 
which involve the extensions of the middle surface, which partly, although not 
entirely, arises from his having omitted the factor 1 + h' ja. It is not improbable 
that these terms may be small, but at the same time we are not at liberty to neglect 
them altogether ; for it is quite evident that a term such as S (2$/x) in the variational 
equation, will give rise to terms in the equations of motion and the equations giving 
the values of the sectional stresses, which do not involve the extension of the middle 
surface. 

In the second place, on comparing Mr. Love's variational equation of motion''^ with 
my equations (25), (26), (28), and (29), it will be seen that he has omitted several 
terms in the expressions for SC, SU, and h%. 

In the third place he states (p. 521) that the extensional quantities '' cti, cr^, tjs may 
not, in general, be regarded as of a higher order of small quantities than /c^, Xj, /c^," 
which are the quantities upon which the bending depends. The argument of Lord 
RAYLEiGHt appears to me to show, that at points whose distance from the edge is large 
in comparison with the thickness, the extensional terms are usually small in comparison 
with the terms upon which the bending depends. It must be obvious to every one, 
that a thin plate of metal or a steel spring can be bent with the greatest ease by means 
of the fingers, whereas the production of any extension of the middle surface which 
would be capable of measurement, would involve considerable muscular efibrL These 
considerations indicate that when a thin shell is vibrating, the change of curvature is 
so greatly in excess of the extension of the middle surface, that notwithstanding the 
smallness of h^ compared with h^ the product h? (8p~^)^ is largej compared with 

^ 'Phil. Trans.,' A., 1888, p. 514, equation (19). 
t ^ Roy. Soc. Proc.,' vol. 45, p. 105. 

J The problem discussed in § 14 shows that the product ha^ may he of the order i^^ (^P~^)^ except in the 
neighbourhood of a free edge ; but in the equations of motion we have to deal with the quantities Jkt and 
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the product ha^. At the same time, inasmuch as the production of change of 
curvature involves some extension or contraction of all but the central layers, and 
consequently of those portions of the shell which are near its external surface, it does 
not seem unreasonable to suppose that in the neighbourhood of a free edge, an exten- 
sion or contraction of the middle surface may take place, which is comparable with the 
change of curvature. 

In the fourth place, Mr. Love appears to have argued as if the equations of 
motion of a shell, whose middle surface undergoes no extension or contraction 
throughout the motion, might be obtained from his general equations (30), (31), (32), 
by putting ctj^ = 0-2 = ^ = ; but it has already been pointed out, that the correct 
equations for this kind of motion must be obtained by taking the variation subject 
to the conditions of inextensibility^ ■ and introducing indeterminate multipliers. It 
will be shown in the next section, that in the case of the flexural vibrations of an 
indefinitely \ot\^ complete cylindrical shell considered by Hoppe and Lord Ray- 
LEIGH,"''" the differential equation for the tangential displacement v is of the sixth 
degree, and that when the cross section of the shell consists of a circular arc, this 
equation contains sufficient constants to enable all the conditions of the problem to be 
satisfied. 

13. The first problem which we shall consider will be that of the flexural vibrations 
of an indefinitely long cylinder, in which the displacement of every element lies in a 
plane perpendicular to the axis of the cylinder, and which has been discussed by 
Hoppe and Lord Rayleigh. 

/In this problem the middle surface is supposed to undergo no extension or contrac- 
tion throughout the motioUp and the solution is most easily obtained by means of the 
general equations (11). In these equations we must omit all the terms on the right 
hand sides which involve A^, for they w^ould, if retained, give rise to a term involving 
/^* in the period equation, which must be rejected, since we do not carry the approxi- 
mation further than h^ in determining the period. 

We evidently havet M^^ z=: N^ = Hg = ; also none of the quantities are functions 
of z. The equations of motion are thus 

^' + N, = 2phav, 

1 m 7 

* ' Theory of Sound,' vol. 1, p. 324 ; ' Roj. Sec. Proc.,' vol. 45, p. 129. Equation (51). 
f We shall presently see that these conditions imply a constraint at infinity. 
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also the condition of inextensibility gives 

Eliminating Ni, T^, and w^ and substituting the value of G^ from the first of (52), 
we obtain 

+ -:z^)v+jT^--v=0 (59), 



3pa^ {m + n) \d(l>^ d^J ^ d(^'^ 



whence, putting 



we obtain 



V = Ae'^^ "^ '"^, 



.3 



^ "^ 3pa^ (m + ^) (s^ + 1) ^ ^' 

which is the required result. 

If the cylinder is complete, s is any integer, unity excluded, but if the cross-section 
of the cylinder consists of a circular arc of length 2aa, s will not be an integer. Its 
values in terms of p are the six roots of (60), but in order to obtain the frequency 
equation, the value of 5 in terms of the dimensions and elastic constants is required. 
The additional equations are obtained from the boundary conditions, which have to be 
satisfied along the straight edges of the shell, and these require that the tension Tg, 
the normal shearing stress N^, and the flexural couple G^ should vanish at the edges 

where ^ = dz ^• 
Since 



where 



^1= ^ 3 (m + ^) ^' 



1 / d^v dv 

^"^ a^ \d4^ "^ d^^ 



the boundary conditions are obviously 

/. = 0, 

dcf> ' 



4-mnh^ d^a . dh 



Spa/^ (m + n) dcjr' d<^ df 

These conditions have to be satisfied at each of the edges of the shell where 
(^ = -J- a, and there are, therefore, six equations of condition ; hence the six constants 

MDOCCXC, — A, 3 N 
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which appear in the solution of (59) can be eliminated, and the resulting determinantal 
equation, combined with (60), will give the frequency."^ 

If a complete cylinder of finite length were vibrating in this manner, it would be 
necessary to satisfy the conditions at the circular ends, and this would require that 
Tj = Qj G^ = at the ends for all values of ^ ; and from the first and fourth of (43) 
we see that this requires that fi = 0, or 

whence 

w = A cos (^ + B sin (jy 

for all values of ^. Since it is impossible to satisfy this condition for the kind of motion 
considered, it follows that when the cylinder is of finite length it would be necessary 
to apply at every point of the circular boundary a tension T^ and a couple Gg of the 
requisite amount. 

This is the question upon which Lord Rayleigh and Mr. Love are at issue ; and 
the preceding investigation shows that Mr. Love is right in supposing that it is 
impossible to satisfy the boundary conditions along the curved edges of a cylindrical 
shell when these edges are free, although he does not appear to have noticed that it 
is possible to satisfy these conditions when the free edges are generating lines. In 
order to obtain a complete mathematical solution of this question, it would be 
necessary to work out the problem of the free vibrations of a complete cylindrical 
shell of given length 21, which is deformed in such a manner that dv/dcj) -^ w = 0, 
where v and w are functions of <^ alone, and is then let go, without assuming that 
the middle surface remains unextended during the subsequent motion. 

Owing unfortunately to the extremely complicated nature of the general 
equations, a rigorous solution of this problem would be exceedingly difiicult. We 
shall, however, be able to throw some light upon this question, by solving and 
discussing the following much simpler statical problem. 

14. Let us consider a heavy cylindrical shell, whose cross section is a semicircle, 
and which is suspended by means of vertical bands attached to its straight edges, so 
that its axis is horizontal ; and let us investigate the state of strain produced by its 

own weight. 

In order to simplify the problem as much as possible, we shall suppose that the 
displacement of every point of the middle surface lies in a plane perpendicular to the 
axis, and we shall afterwards .investigate the stresses which must be appHed to the 
circular edges, in order to maintain this state of things. 

* [Tkis problem is of a similar character to that of a bar, whose natural form is circular, and which 
has been discussed by Lamb. * London Math. Soc. Proc.,* vol. 19, p. 365, — June, 1890.] 
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We have 



Y = — ^ sin ^, Z = g cos <j& ; 



whence, if W = 2gpah, the equations of eqiuHbrium are 



from which we obtain 



the integral of which is 



d(j> 
d(j> 
d(l> 



+ Ni = W sin ^ , 



ig — 



W cos ^ , 



h^ 



4- N^a == -— r-- W sin ^ ; 



?>a 



cPT 



3 



d^ 



3 



+ T^ = 2W cos (f> , 



and, therefore, 



Tg = A cos ^ + B sin ^ + W^ sin </> , 
Ni = A sin <j& — B cos — W<j& cos ^. 



Since N^ = when ^ = |^ tt, A = ; also since T^ = | ttW when <j& = 
whence 

and, therefore, 

dG^ 



i^^ 



W^ sin <j&, Ni = — W^. cos ^ 



(i0 



^2 

Wa {cl> cos (^ - ^2 sin </>), 



3a' 



whence 



^2 



G^ = Wa (<j& sin ^ + cos </> + — cos (j)) + 0. 



3^2 



(61), 



Since Gi = when <f) = ^ tt, G = — ^Wira ; accordingly 



Gi = Wa ^ 



h^ 



^ sin .^ + 1 + — COS <j& — 4ir 



3a^. 



3 N 2 



(62). 
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G^= -in¥:^= -inh^l + E)|Ji . 



ft • • 



^nh^lJL 



Wa 
1 +E 



<j) sin ^ + ( 1 + ^ ) cos 6 



2 



IT 



• * e 



■ (63), 



(64). 



Again, if R denote the change of curvature along a circular section, so that 



R= - 



1 /^ 

a2 \ d(j)^ 



+ w 



we have 



Also by (18) 



ju, = II — Ha-^/ct 



9 9 



, (65). 



/^ 



- (2 + E) ja/a + ~~^ 



E c^2 



OTc 



a^ d(f>' 



and, therefore, 



d?'=z=(l + E)iic'= -(1 + E)(2 +E)/x/a + 



(1 + E)E 6^V, 



3 



(X-^ 



(i^' 



whence, by the second of (44), 



T, = inh (1 + E)o-, - '^(1 + E)(3 + 2E)/. + 5^ {1+ E) E^ + ^.EWcos,^. 



Substituting the values of T3 and /x from (61) and (64), we obtain 



2«^ (1 + E) 



20-3 + 



hmd? 



(To 



Za? d^\ 



+ 1(2 + 3E) ^ ^ sin (^ + M + ^ ) cos <ji — l tt 

+ --^ W cos (^ — W^ sin ^ == 0. 



This equation might, if necessary, be solved by successive approximation, but a 
first approximation will be sufficient. Omitting the terms in h?, and recollecting 
that W involves A as a factor, we obtain 

4:nh{l +E)crg+ W(cos(^ — Itt) + | EW (</> sin <^ + cos <^ -»• ^ tt) c= . (66), 

whence from (64), (65), and (66), we obtain 



E E 



a 



+ 



'da (^TT — (j> sin <f) — cos <jb) 
/r { J TT — cos <5b + f K (^ TT — ^ sin ^ — cos «^)} 



* o 



(67). 
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Since the numerator of this fraction is an even function of ^, it does not change 
sign with ^ ; also the numerator is always positive between the limits -| tt and — |- tt, 
and its maximum value occurs when ^ = and is equal to ^tt — 1, and its minimum 
value occurs when ^ = |^ tt and is equal to zero. We, therefore, see that when ^ = 0, 



and when (^ = |^ tt. 



0"3 


E 


E 

a 


+ 
E 


Sa 




<^2 




a 


» 



Since the thickness of the shell is supposed to be small compared with its radius, it 
follows that the change of curvature is large compared with the extension of the 
middle surface, except when a (Itt •— ^) is comparable with h, i.e., in the neighbour- 
hood of the straight edges of the shell ; and therefore at all points of the shell whose 
distances from the edges are large in comparison with its thickness, the terms 
depending upon the product of the change of curvature and the cube of the thickness, 
^.e., the terms upon which the bending depends, are of the same order as the terms 
depending upon the product of the extension of the middle surface and the thickness ; 
but at points whose distances from the edge are comparable with the thickness of the 
shell, the extension of the middle surface is of the same order as the change of 
curvature, and therefore the terms depending upon the product of the change of 
curvature and the cube of the thickness are small in comparison with the terms 
depending upon the product of the extension and the thickness. 

We shall now calculate the stresses which must be applied to the circular edges in 
order to maintain this particular kind of strain. From (43) we have 

Substituting the values of ju, and cr^ from (62), (63), and (64), we see that the terms 
in <r^ may be omitted, and we obtain 

Gg = -z IT, (1^ TT — ■ <j& sin ^ — cos ^) ..... . (S8), 



1 + E 



which shows that G^ is positive. 
Also 



Ti = 4^^ Ecr, --. ^' E(2 + 3E) /x 
EW 

which shows that T^ is positive. 

Comparing (68) and (69) with (61) and (62) we see that ratios of the tension 
Tj and the couple Gg, to T^ and G^ are numerically equal to E/(l + E) ; we further 
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see that G^ is negative, and, therefore, the strain tends to increase the curvature of 
the circular sections. Now when a cylindrical shell is bent about a generating line in 
such a manner that its curvature is increased, all lines parallel to the axis which lie 
on the convex side of the middle surface will be contracted, whilst all such lines 
which lie on the concave side will be extended, and this contraction and extension 
will give rise to a couple about the circular sections which tends to produce anticlastic 
curvature of the generating lines. In order to prevent this taking place it is necessary 
to apply at every point of the circular edges a couple G^ tending to produce synclastic 
curvature, and a tension T^ parallel to the axis, whose values are given by (68) and 
(69). If this couple and tension were removed, the middle surface would bend about 
its circular sections, and anticlastic curvature of the generating lines would be pro- 
duced, and this would necessarily involve extension or contraction parallel to the axis, 
so that the problem could no longer be treated as one of two dimensions. 

It must, however, be within the experience of everyone that when a thin cylindrical 
shell of finite length, whose cross section is the arc of a circle, is bent about its 
generating lines, the shell does not assume a saddle-back form, and consequently the 
anticlastic curvature of the generating lines must be so small as to be inappreciable. 
This circumstance furnishes an additional argument in favour of the supposition that 
the extension of the middle surface is only sensible in the neighbourhood of the free 
edges. 

We therefore conclude that if the circular edges were free, some extension or con- 
traction of the middle surface must necessarily take place, but that this extension or 
contraction is small compared with the change of curvature along a circular section, 
except just in the neighbourhood of the edges. From these considerations the infer- 
ence Is, that if by means of proper constraints applied to the circular edges, a 
cylindrical shell were enabled to execute the non-extensional vibrations discussed in 
§13, the vibrations would cease to be non-extensional if the constraints were 
removed ; but that the amplitudes of those portions of the displacements upon which 
the extension depends, would be very small compared with the amplitudes of those 
portions upon which the change of curvature along a circular section depends, except 
just in the neighbourhood of the edges. Moreover, the theory of plane plates shows, 
that the frequency of the extensional vibrations is expressible^ by means of a series of 
even powers of h, commencing with a term independent of h ; whilst the frequency of 
the flexural vibrations is expressible by means of a similar series com^mencing with h^. 
It therefore follows, that the pitch of the notes arising from the former class of vibra- 
tions, is high compared with the pitch of those arising from the latter class. And 
although, except under special circumstances, it is not possible in the case of curved 
shells whose edges are free, for these two classes of vibrations to coexist independently, 
as in the case of a plane plate ; yet recent investigations show, that the pitch of 

* Lord Ratlbigh, ' London Math. Soc. Proc.,' vol. 20, p. 225. See especially equations (38), (45), 
and (53). 
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notes which mainly depend upon the extension is usually high in comparison with 
the pitch of notes which mainly depend upon bending, and consequently the 
notes arising from the former cause, both on account of the smallness of their 
amplitudes and the highness of their pitch, would probably be so feeble in comparison 
with those which arise from the latter cause, as to be scarcely capable of producing 
any appreciable effect upon the ear. Judging from the usual course of such 
investigations, the probable form of the exact solution of the problems suggested 
at the end of § 1 3 would be that of an infinite series, the periods of the different 
components of which would satisfy a transcendental equation having an infinite 
number of roots ; but the preceding considerations point to the conclusion that the 
frequency of the gravest"^ note given by (60), viz., p'^ = 4:8mnh^ /5 pa"^ (m + n), 
although perhaps not rigorously accurate, is a close approximation to the truth. 

Spherical Shells, 

15. The fundamental equations for a spherical shell can be investigated in precisely 
the same manner as in the case of a cylindrical shell. 

If u\ v\ w' be the component displacements at any point of the substance of the 
shell in the directions, ^, ^, r, the equations connecting the displacements and strains 
are 




CTi = 



1 Idu' 
't\M 



+ w\ 



1 dv^ 



or 



3 



(To = 



T \sin 6 d<p 

dvJ 
dr 



+ u cot 6 + '^0 J 



TJS 



m 



1 — 



2 — 



1 dw^ dv' 

— - ■ — 1 1 1 ' 

?' sin dcj) dr 

d%b' u' 1 dw' 
dr V T dO 

1 fdv' 



V' 

r 



( ' ' ' • • ' X^/i 



TXT 



= :At^ -t;^cot^+ t 



1 dtt^ \ 



dO 



sin 6 d(f> 



* [The experiments of Lord Raylbigh, * Phil. Mag./ Jan.j 1890, show that the effective pitch of a 
bell is usually not the same as that of its gravest tone ; and, in the bells which he examined, the fifth 
tone in order was the one which agreed with the nominal pitch of the bell. — June, 1890.] 
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whence 



also 



^dr^ ^ 

'd^w" 
dr' 



dr 

'dw 
dr 

'dw' 
dr 



^ 4-- 



^ a 
A 



m + 7h 



2 

dr I 
dr 

m + n 



2\ I "^3 



a 






1 dvj 
a dd 



~\ 



dw 



a sin 6 defy 



y 



EK 



^ 



dA E dK 

' ' '*' . — .... ■ — . - M,| L,Mi — .— __ 

a (m + 7i) dO a dO 



' ■ (2j, 



"^ 



, vsy 1 <^-A , E dlL , 

a a {m + ?^) sin ^ d(^ amxQ d<^ f 



* » \^)* 



-- E (X + ja) 



16. We can now obtain the equations of motion in terms of the sectional stresses. 

If dB be an element of the middle surface whose coordinates are {a, 0, (f)), and 
dS' an element of a layer of the shell whose coordinates are (a + h\ 0, (})), then 
c^S' == (1 + h'/ay dB ; whence, if m the figure OA, OB respectively coincide with the 
meridians and circular sections, we obtain by resolving parallel to OA, 



"^ (T^a sin e S(^) S^ - T^a cos 80 Scjy + |^ (M^a 8^) Sc^ + N,a sin 80 S(^ 



dd 



d^ 



= p dS {' fu' -- X) (I -\- h' /af dh' . . (4). 



But 






dr 



accordingly if we substitute the values of {du/dr) and {d^u/dr^) from (2) and (3), and 
recollect that all quantities which vanish with h may be omitted when multiplied by 
h^^ the right hand side of (4) becomes 



h^' 



pdB \2h[l + ~2 P^ + 



a' 



da d0 "^ 3a^ d0 



2h(l + 



h? 



?fa' 



X 



Kesolving parallel to OB, OC, and then taking moments about OA, OB, 00 we 
shall obtain in a similar way five other equations, which, together with (4), may be 
written. 
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d 



dM 

(Ti sin 6) - T3 cos ^ + -^ + N3 sin 



■^ 



2h 



!+!)» + 



d(p 
hm dK 






(X-* 



3a dO U^ dd 



— 2^^ fl -{-. r-^ j X >pa sin (9, 



■54-.4(M,.in.) + M.o„.. + K..„. 



d0 






4Jh^ dw 



a" 



oct sin 9 d(j> da^ sin 6 d^ 



2h(l + 



3a^ 



Y >pa sin 6^ 



da 



^ (N, sin ^) + ^ - {Ti + T,) sin ^ 

2/^/1 + ^A ^(, _ -L AS E (X + ^) _ *;f^ EK 

"7~ + Nja sin 6 + ^^ (Hi ™i ^) "" H^ cos 6 
= f ph^ Iz^. 5 - 8v + 2Y^ sin^, 



-- 2h(l + Tl) ^ f P* sin ^, 



I. 



sin ^ d6 



d 



(Gg sin ^) + Gi cos ^ — NgCt sin ^ + 



C^H: 



= _ I p/,3 /|^ _ 3tt + 2X) sin 0, 



(Ml — Mg) a — Hi — .Hg = 0, 



(5) 



17. We shall now (as in the case of a cylindrical shell) proceed to obtain the 
values of the couples and the stresses M^, Mg by direct calculation. 
We have 



whence 



k 



Tj a sin 6 S(j> = P' (a + h') sin 6 8^ dh' ; 



~/« 



T: 



T,= 






6^r^ 



3a \ dr 



Ml = M^ = 2n;itir + i n¥ 



'd^ 



GTc 



(ir^ 



2nl\ Idxss^ 



Gi 



Gg = 



3 '^ 



dr I a 



2.^3 
3 ' ^ 



dr 



+ 



p 



a 



Hi 



TT 



%nU 



3 



'^ 



ttTf 



Wf 



MDCCCX.C. — A. 



dr 
3 



1 






"1 



> 



. . . (6). 



J 
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The third and sixth of these equations satisfy the last of (5) as ought to be the case. 
Also, employing our previous notation, we see that 



G,=:-|r./.8 d? + 



a 



i, C^^ = %nhP 



H,= -H,= -|n;.3(^^ + ^3^ 



+ 



a 



> 



■ ■ in 



_>' 



Since the couples are proportional to the cube of the thickness, it follows from the 
fourth and fifth of (5) that the normal shearing stresses N^, Ng are also proportional 
to the cube of the thickness, and, therefore, that the shearing strains tjt^, zj'^ are 
quadratic functions of h and h\ 

Employing our previous notation, the next thing is to calculate the quantities 
X, /x, p, \\ II , p. We have 



X=: 



/x 



dr 



1 fd^w 
a^Vdf^ 



E 



+ ^) --(^1+ 0-2) 



'd 



CTc 



P 



dr 

'dnjr^ 
dr 



d^w ^ dw , 

a^ \sin'^ 6 d(p^ dO 

^ dv) d^vj 

(^ sin d\ d(f> d6 d(j)^ 



E 



a 



~] 



{<r, + <T,) y . . (8), 



in which equations we have omitted all quantities which vanish with A, because 
X, fi, p occur in expressions which are multiplied by h^. Similarly 



X'=-^- 



^ = — 



p = 



a 
a 

a 



E , ' X . E d^K 



■^ 



E 



a 



(X + /x) + 



^3 d6^ 
E/ 1 



a^ \sin^ 6 d(j)^ 



d^K ^ , . dK 
+ cot u 



d0 



r 



. (9). 



2E 



a^ sin 



, ^dK #K 
cot u -r: — 



d(f) dO dj) 



18. The variational equation may be written 



8W + 



SU + SIL 



(10), 



and we must now calculate the values of the four terms in it, and we shall begin 
with W. 

Since w^e may omit ct'i, Tn\, and may, therefore, write t?y for ^3, the potential energy 
of any portion of the shell is 

W =^ i fff' [(m + n) A'^ + n {^'^ - 4. (cr\a-',+ cr>'3 + o-Vi)}] (1 + h'/aydh'dS (11) 
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where the integration with respect to S extends over the middle surface of the portion 
considered. Since 



^' = ^ + *'(f) + i*'^(S)+ • ■ ■ 



we obtain 



also 



and 



h 



|(m + fO ^''(1 ^h'/afdh' 



h 



— — \lflh *~j" Th) 



M'+aA^+iWfT+iW^)+f-(f 



m -\- n 



j ^ (l + ^) (o-i + 0-2 + A/2n)^ + -| W (X + /x) 



2 



+ i ^^ (o-i + o-g) (X' + f^O + -^ (o"i + ^2) (^ + /^) K . . . . . (1-2), 



2n 



h 



Ji 



(T\a'\ (1 + h'jaf d¥ = 4nh (1 + -| h^/a^) c^o-c^ + | n/i^x^ 



8nM 



+ I nA^ (XVg + fi'o'^) + -^^ (Xo-^ + ixo-^) . . (13), 



2n 



k 



^k 



{a-'i + cr'2) o-'s (1 + KjafdJh 



= 4:nh{l + ihya^) 



A 



m + n 



E (cr, + a,)] {a, + cr,)^ f ^^^1 (X + i^f 



f t^/^^E (o-j + 0-2) (X' + ii) 



l&nM 



E (X + /a) ((T]^ 4" ^%} (14) ; 



lastly 



2 



n 



f w'^ (1 + h'laf dN =znh(l+l h^a^) w^ + i nhy + ^n¥mp + 



--z — mp (15). 
da -^ ^ ^ 



Substituting from (12), (13), (14) and (15) in (11), the value of W per unit of area 
of the middle surface is, 

W = 2nh (1 + ^) {a,' + cr/ + E (cr, + a.f + i ^'} 
+ t n¥ {X2 + ^3 ^ E (X + /.)2 + I _p3} 






3a 



3 o 2 



• * » 



» ♦ • » • (Xu/, 
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We must now obtain 8C. We have 



S€ = p\\W (u' Sw' H- v' Sv' + io' Stv') (1 + h'/af dh' dS ; 



al 



so. 



'k 



-h 



h^\ '* 



du dSio 



d^Sto . dH 



u^ Su' {I + h^af dh^ = 2hll +--~,)uSu + ih^ ^ 



3d 



dr^ dr^ 



, AW' [•■ dBu , du „ 
3a \ dr dr 



-, , ]i^\-- ^ . 2¥ fdw ••\[dhw (, 



tiay 



3a^ \d$ 



\de 



hmf--dBK , dK^ \ AW 



Ja \^ dO '^ dd 



Za' 



dSw ^ \ , [dw . ^ 



d0 



dO 



by (2) and (3). Treating the other terms in a similar way, we obtain 



8C = 2p/i (1 + i h^la^) ' (u 8u + v8v + w Bw) dS 



+ IM' 



1 /dto ' '\ fdho 
a^\d0^^^J\dJ 



Sii) 4 



r 



• * 

1/1 dw 

^ \sin 6 dcf) 



3v 



1 dSta 



sin 6 d^ 



Sv 



+ EfEK-^^^^SKJc^S 



+ ipIM 



r r 



U 



dSK 



+ 



V 



dSK 



1 dK. 



\a dO a sin dcf} 



iv{S\ + Six) + -^jSu + 



1 dK 



a sin 6 d^ 



8v 



• • « • I 

(X + ft) Sw > dS 



4cpW 
3<x . 



1 (^ _ u) 8u + U ^, f - v) Sv + EK 8w 

a \dd / a \sin U d<p j 



d% 



* • 



. . (17). 



We must next find S%, 
We have 



811 



'k 



■h 



(P' Su + W 8v') {a + h') sin 9 dK d4> + 

+ N^a sin 68wd(l> + 



'k 



h 



' Si;' + W hu) {a + ^0 dJi dd 



N^a 8iu dd 



» • 



(^loj. 



whence 
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T, Su + M, 8v + N, Stv 



Gr.3 /dSlV 



8'u] + 



Hi / 1 d^o 

a \sin 6 dcf) 



Sv 



. 27ihm^dSK , nhmuj dSK] . . , , 
^ 3a dd ^ Sa sin 6 d^ } ^ 



+ \\M,8u+ T^Sv + N, Sto + ^ r^ "^ - 8^ 



Ho fdSia 



'3 






<x \siii (^(/) 



a \d6 



m 



Lastly, 



8U 



Ch 



,h 



^ Sa sin (9 d^ ^ 3a ^6> J 



(X 8u + Y Si;' + Z 8^(^0 (1 + Ii/af dh'dS 



« • 



2phll+ £-) If (X 8tt + Y S« + Z Sit;) (iS 



+ ip/^'E 



~3a 



X dSK Y (iSK r7 /C>x I Q' \ 

[a ~dJ + ^^iiT^ "rf^ ~ Z {8X + 8/^) 



dS 



X /(^S^^; 
a \ dO 



^-)+!(sib'^-^^) + ^^^K 



(19), 



• d^ . . (20). 



19. We shall,as in the case of a cylindrical shell, denote the four lines of W by W^, 
W3, W3, W4. Whence 

SWi = Anh. ( 1 -h i h^la^) f | (^So-i + ^^o-^ + k ^^^) ^^ 



inh ( 1 + 



h^ 



\ ■ Za^ 



(^8i( -\-\vslv)a sin B d^ -\- (^ hv -\- \ ns Sm) a d^ 



- ^^^ i^ + Ir^) 



Ad 



(^ sin ^) - 2$ cos ^ + i ^ j Stt 



(l^\ 



+ -! ^+ i Ta (^ sin ^) + 1 OT cos e \lv - (ft + 3$) sin ^ 8w 



(^^ 



c^d" 



addd(f) (21); 



from which we obtain the approxima,te equations 



Ti = 4«/iE, T3 = 4«^2S 



;^ = Mg = 2n^KT 



(22), 



p'« = 



pv = 



2??, r d 



asm 6 \dd 

2n (dM . , d 



(^ sin ^) - B cos 6' + 1 



dm 
d(j)^ 



-f- pX. 



a sin ^ \ (i^ 

271 



+ l^(t;ysm^) + -|T^cos^l +/>Y y. . . (23) 



p^^ = (^ + 3$) + pZ 



470 



MR. A. B. BASSET OK THE EXTENSION AND FLEXURE OF 



These axe the equations which have been obtained by Mr. Loye/' and which have 
been employed by him in discussing the extensional vibrations of a spherical shell. 
Again 



SW^ = f nP 



Substituting the values of X, /x, p from (8) we obtain 



|[3ESX(^S = 



HMBlHaM 



also 



de^ "*" '^ \de • mud d<j> 

d 



(dhh , 1 dhv , ^v , /, . ^cs \ . /I 7/1 7 , 

— + zrrix ~7T + o'?^ cot a + 2^^t; ) ^ sm (^ do' rfc^ 



:"?Si^ 



E|E sin ^St^ — ^ (iE sin ^) S^(; + 12 sin ^ ~j~ I dcj) 



d6 

' ^ . d% d^ 

EsinO-rz Su 4- 'E-r-Sv 

dd d^ 



de 



ElESt' dd 



d} 



de^ 



(igsin^) + {l + 2E)iEsin^lS 



w 



d0d^ (24), 



f|#S/xrfS = ^ jjdp 



1 d^Sw . ^ dSw , ^ . ^ 

+ cos c7 ~j- + dt{; sm ^ 



sin d(^ 



3 



c^(9 



+ E ( -j^ sin ^ + ~jj -f S^^ cos ^ + 2S^^ sin (9 ) !> (i<9 (i<j& 



de 



dcf) 



Wn - niw — 



(E;fp sin ^St^ + d?" cos ^Stf;) d(f> 

sm ^ (^^ sm 6 d^ ) 

Esin^^St^ + E^St; 



1 d^§ d 



sin d<^' 



de 



(# cos + (1 + 2E)#sin e \ hiv 



ddd(f> (25). 



In the last termp Sp, we must treat the integral which involves d^Zwjdd d(j) exactly 
in the same way as in the corresponding case of a cylindrical shell, and we shall thus 
obtain 



2 



pSpdB=: 



I , ^ dZw 



m) '" '^ 



1 [i^ 



dS 



w 



dp , d^p 



cot^ 



d<f) 

ammmmmm 



< 



p cot ^ + 2 m ) ^^ *"" iP 



dS 



w 



dej 



dd 



d(j> d6d^^ 



Swd6d(f> ......... (26). 



* *Pliil. Trans.,' A, 1888, p. 527. Equation (23) coi*responds to Love's equations (46), (47), and (48) 
and (22) to (72). 
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Adding (24), (25), and (26), we finally obtain 



sw 



2 



~3^ 



- {iE + E (iE + dF)} 8u -ipSu+ -; 



d 



sin 6 dd 



+ s 



, dp 

X - 

dj) 



> Sw 



'dSiv ^ \ . /I 



(m sin 61) - d?" cos 

dSw 



sin ^ ^<^ 



— Sv 



a sin ^ defy 






ii>S7^-.{d?+E(iE+dr)}s^+-: 



-# 



sin \d<p 
1 (iSi^ 



^ +j9cos^ + |-sin ^ -:^] St(; 



de 



^sin c?(^ 



S^j - i P (-^ - S^^M ^ ^^ 



+ -f ?l/?y^ 



E sin ^ I (5£ + dF) S« + E ^ (S + dF) 8« 



-{|^(i£sin^)+(14-2E){IS + #)Bin^ + ^^^^-|(dFcos^) 



+ 2(l+E)dFsin^ + cot4^ + J^^S^ 



c?^c?<^ (27). 



The expressions for Wg, W^ niay, as in the case of a cyhndrical shell, be divided into 
two parts Wg', W3", W/, W"/'. The values of SW^\ SW/, may at once be written 
down from (21) by changing ^, IS, ttt into IE', dF'? p' ^^nd iS, dF? P respectively, and 
by altering the coefficient into -| nh^ and 8nh^/Sa respectively. With regard to W3'' 
we have 

SW3" = I nh^ f f (^ SV + 23 S/x^ + I TIT 8p) dS. 
Substituting the value of V from (9) and integrating once by parts, we shall obtain 



If^ SX' c^S = E f ^ ^sin e # 



ti ti 



dO 
E -^ (^ sin ^) "^ + fta sin ^ { 2 SX + E (SX + S/x) } 



dO 



dO 



dd d(j>. 



Treating the other terms in a similar manner, we shall finally obtain 



SW3" 



2rfE^6ZSK nhm^ dSK 

+ 



3f^ 



d0 3a sin d<f> 



a sin d(l> 



+ 



2^ASEa3^SK nhmjjjdSK 

+ 



36t sin d^ 



da d0 



add 



+ 



27iM 

E 1. 



E 



sin d0 



'^ (^sin^)-2Scos^+l^5"' ^^'^ 



sm' 



H#+*S'^^^+""^'^ 






d^\ d0 

+ Ea (^ + 3$) (SX + S/i) 



c^S 



- ^^f{(^8^ + 3^S/^ + it=v8p)rfS. 
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If in the first surface integral in this equation we substitute the approximate values 
of the coefficients of d hlLjdd &c. from (23)3 which we may do since this integral is 
multiplied by ffi, and. then substitute the vakies of 8W3'', SC, h%, and SU in (10), it 
will be found that all the terms involving d S^K/dO, d SK/dcj), and SX + S/x cut out ; 
we are therefore no longer concerned with them, and the value of SWg'^ reduces to 
the last line ; on this understanding we may write 



SW3" + 8W/ =: 



3 a 



(^ 8X + ^ 8fc + I to- Sp) cZS 



. (28). 



The variation of the right-hand side of (28) might at once be written down from 
(27) by substituting ^, M, and trr for IE, ^, and p ; but it will be more convenient 
to present the results in another form. Taking the first term, and integrating the 
second differential coefficients once by parts, Y7e obtain 



8X6^8 = 






sin d(j) -^ 



~,{^sine)'''' 



dd 



de 



sm 



dZ 



■w 



- E^cKsin^8E| ded<l>. 



Treating the other terms in a similar way, and adding to the result from (21) that 
portion of SW^^ which depends upon h^, and finally replacing the coefficients of 
d Sw/dO — Su, &c., by their approximate values from (23), the final result will be 



4:nh^ f f ^ fdStu 



\ 



Su) + 



2 



tjy 



1 dSw 
^sin 9 d(j) 



4znW 



\de 

I 1 dhw Q. \ , 1 (dZw 

(u - X) (^1^ _ 8«) + (;; - Y) 



— hv\\ a sin 6d^ 



hit] \ add 



1 dhw 



^sin 6 dip 



Sv 



+ Ea (iv -- Z)8K\dS (29) 



This result enables us to test the accuracy of a portion of our work, and the funda- 
mental hypothesis on which the theory is based; for if we substitute in (10) the 
expression (29), and also the value of 8% from (19), it w^ill be seen that we have 
reproduced the values of the couples which are given by (7) ; also comparing with SU,, 
the line integral parts of SW^, given by (27), the fine integral parts of SWg^and SW/, 
which, as we have explained above, are obtained from (21) by changing certain letters, 
we see that we have also reproduced the values of M;,, M^, given by the third of (6). 
We may, therefore, omit the couple terms, and also the terms in M ; also, since we 



CZLIFDEIOAL AFD SPHEEICAL THIN ELASTIC SHELLS. 



473 



have disposed of the terms in SW^, which involve h?, we shall write SW^' for the 
remaining portion which depends upon h, and the variational equation finally 
becomes 



^W( + SW/ +. SW3' + SW/ + 2p/^I + i Wla^) {;u hti + vhv + w hw) dB 

• • • • 

1 dw _• A / 1 dSw 



+ fM 



3 



1 /dw r'\ fdSw cs \ , 1 
-- 2u I { —7:; bu]+ -I 



a^ \dd 



dd 



a^ \sin 6 dj> 



2v 



Bin 6 d<f) 



Sv 



w 



+ E(EK--^)SK 



+ iM«E 



ApW 



IdK 



d^ 



• • • • 



%/ J 



- — r- S% H r 7, 77 hv — (X + u) S(^ 

a «xp a sm d d(f> ^ ' ' 



dQ 



'Sa 






c^S 



2ph 1 1 + 



It 



X 



XSii -{- Y Sv + Z Sw) dS 

1 dSw 



^dSw 
'dd 



8u) -{-Y 



,sin e # 



-bv] + ZEaBK 






dS 



+ (Tj Sw + N3 §w) a sin 6 d(l) -\- (T^ 8t; + Nj 8iv) a dd 



dS 



« • • « 



(30). 



We have now got rid of all the terms involving the second differential coefficients 
of Su, Bv, 8w ; and all that remains to be done is to integrate by parts the terms which 
bvolva the firrt differential coefficient.. Putting 



a 



dw 



in + X, /8 



1 dw 



sind d^ 



- 2v + Y, 7 = E (aEK -w-\-Z) (31), 



we have 



2p¥ 



ff 



dSta , ^ dhw , c^TT^ Id 



dd Bind d4> 



2p¥ 
3a , 



(yStt + f^Bw) a sin ^c?^ + J^\ ^y§^ ^ /3S^) adl^ 






• (ysin^) — ycos 0>Zu-\-~~Zv'\'\ jz{ami ^) + — — 2'}/sin^ }►S^c^ 



de 



dS 
sin 



(32), 



Substituting the values of SW/, SW/, SW/, SW/, and the right hand side of (32) 
in (30), and picking out the line integral terms, we obtain the following equations 
for the sectional stresses, viz., 

MDCCCXC, — A, 3 V 
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Anh^ 



1\ = inhWL + -l^m -.E(|g 4-#)} -hinhm' +^E(aEK - w + z) 



2 — 



2nhnT + ^^r P + 3 '^^h^p 



N3 

Hi 



4:nh^ 



3 a 

d 
dd 



^(iHsin^)-.d?cos^ + i^Sl + ^Y^^>^2;i+X 



3a sin 6 
I n¥ (IE + ^/a) 
— f nA.^ (p + Ts-/a) 



d^ 



3a \d9 



> (33) 



which give the values of the sectional stresses across a parallel of latitude ; and 



M;l = 2nhrs + -;^7- p + i w/i^p' 



-N 



T^ = 4.nhU + 



3a 



3a 



{dF- E (IE + dr)} + in¥:^' + ^'E(aEK-w + Z) 






4«^^ ''^^ + l>cos^+isin^^)+^(;^,~-2« + Y 



3a sin 0\d(j) 

- I «/i'^ (# + 3S/a) 






2p7z.^ / 1 dio 
3a \sin ^ d(j> 



Hs = I ^/i^^ (p + trr/a) 



^ (34) 



J 



which give the values of the sectional stresses across a meridian. 

In these equations we may, if we please, substitute the approximate values of 
u, V, w from (23), and by means of these values it can be shown that the values of 
Ni, Ng agree with the values which are obtained by substituting the values of the 
couples in the fourth and fifth of (5). 

Picking out the coefficients of Su, Sv, Sw, in the surface integrals, we obtain the 
equations of motion, which are 



2 ( 1 + :i ) ^^^ + 



AJi^ dw 



V- \ 



a- 



3a dO 



3a^ dO 



^(^+£>^ 



pa sin 6 



=: An < 



+ 



dd 

~3a} 
3a 



^(^sin^)-3$cos^ + l^ 



E sin 4 (IS + dF) + f nh^ 1^ (IE' sin 6) - dP' cos + J- 



dO 



dd 



d^\ 
dcj^j 



d^ 
dd 



(iEsin^)^#COS^ + i 
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= An \ —rr + i ;:^ (^ sin 6) -\- ^zs cos 6 



d^ 



dO 



Anh^ E 4 (IS + d?) + I n¥\'^ ^-h^AP sin e)->r\p cos 6 



+ 



8a 



So- 



(i^ 



c:^*^ 



de 



L# + ^ ^ ^^ ^'^ ^) + ip cos ^ 



2ph? ( 

0<X \ 



+ ^ { /Q sin ^ + 






(36). 



2(1+ —^ (w - Z) - i^^E (X + iti) - ^ EK 



3a7 



3a. 



joa sin 6 



= - 4« (fl + J$) sin ^ + ^ j-^^ (52 sin ^) + (1 + 2E) ( 



+ :i 



1 d^§ 



sin ^ <^^^ 



(^612 



A (j^ COS m + cot ^ ^ + yf-^" 
d6 ^ ' d(}> ' d$ dj) 



+ .-fF) sin ^ 



- 1 wF(iE' + dF')sin d - ^~ {m + #) sin e 



+ 



3a 



3a 



^{asin^)+^-2yBin^ 



(37). 



The correctness of these equations may be tested by substituting the values of the 
sectional stresses from (33) and (34) in the first three of (5), when it will be found 
that we shall reproduce (35), (36), and (37). 

.20, The boundary conditions for a spherical shell may be investigated in exactly 
the same manner as in the case of a cylindrical shell, by means of Stokes^ theorem ; 
for in the present case the theorem may be written 



( TT oW + H. 



d<p 



d(j> 



dc^ + \(^^Sw^K^^]d9=0, 



de 



dd 



the integration extending round any curvihnear rectangle bounded by two meridians 
and two parallels of latitude. If, therefore, in the figure we apply to the side AD 
the stresses, 



to the side BD the stresses 



a sin 6^ d^ ' 



M/ = H7a, 



N/ 



1 dE! 

a de ' 



Hf 
. 2 



= -W) 



and to the sides OB, OA, corresponding and opposite stresses respectively, the 
preceding integral becomes 

3 p 2 
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M/ Bv + N; 8tv + -'- 1-.^-". f J' - 8v)\a sin d4> 



+ \ M/ Zu + Ni' 8t£^ 



Ha'' fdZw 



d9 



— Stn I acZ^ = 



a 



which shows that the work done by this system of stresses is zero. 

If, therefore, we suppose that the rectangle OADB, instead of being under the 
action of the remainder of the shell, is isolated, and that its state of strain is 
maintained by stresses applied to its edges, then it follows that if instead of the 
torsional couples H^, H^, due to the action of contiguous portions of the shell, we 
apply torsional couples ^i, ^3, where 



Hi4- 



(38), 



ra 



H„-H' 



(39), 



the state of strain will remain unchanged, provided we apply in addition the stresses 



1 dM! 



% 






>• 



a sin 6 d^ 



* r 



. , . . (40) 



a 



nd 



1 dUL 



Jti = N, + 



> 



a dd 



• • » • • . ( 4t A f. 



whence eliminating H' between (38) and (40), and between (39) and (41) respectively, 
w^e obtain 



'2' 



a sin 



d 



fi = Noa sin 6 

dji ^ 



and 



d0 



d<j) 



> 



. . . . (42) 






^1^2 _ XT ., , ^1^2 > 



dd 



• • » • « 



. . (43) 



in which we are to remember that |^i =: — l^^, and Hjl = — Hg. 

In these equations the Roman letters denote the stresses due to the action of 
contiguous portions of the shell, whilst the Old English letters denote the values of 
the actual stresses applied to the boundary. If, therefore, the shell consists of a 
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portion of a sphere bounded by two meridians and two parallels of latitude, and 
whose edges are free, the boundary conditions along a parallel of latitude are obtained 
by equating the right hand sides of the first and fourth of (33) and of (42) to zero ; 
whilst the boundary conditions along a meridian are similarly obtained by equating 
the right hand sides of the first and fourth of (34) and of (43) to zero. 

21. If the shell is supposed to vibrate in such a manner, that its middle surface 
does not experience any extension or contraction throughout the motion, the equations 
of motion can be obtained by taking the variation subject to the conditions of 
inextensibility, and introducing indeterminate multipliers. 

22. It will now be convenient to make a short digression for the purpose of con- 
sidering some of the quantities involved. 

Let P be any point on the deformed middle surface whose undisplaced coordinates 
are (a, 6, ^). The coordinates of P after deformation are 

Ii=a + w, @ = 6 + u/a, <l> = ^ + v/a sin ^ . . . (44), 

and since u, v, w are functions of 6 and ^, the elimination of the latter quantities from 
(44) will give a relation between R, @, #, which is the equation of the deformed 
middle surface. 

If pi be the radius of curvature at any point of a meridian section after deformation, 
and P the perpendicular from the centre on to the tangent at that point to the 
deformed section, 

1 — i^ 

Now 



W 



1 + 



(dwld0) 



2 



W{l + du/add)\ 



>i 



and therefore, neglecting cubes of displacements. 



Also 



1 /dw\^ 



dd 



whence 

11 1 fd?w 



Pi 



a~ a^ \de^ 



-{-w] ....... (45). 



which gives the change of curvature along a meridian. 
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We shall now find an expression for the change of curvature along any great circle 
which makes an angle y with a meridian. 

In the figure on p. 463 join OD, and let the angle OED = x? ^^^ ^^^ angle 
DOA = y ; then by (45) the change of curvature along OD is 

1 /dhv , \ 

If w + Sw be the normal displacement at D, it follows by equating the two values 
of Sw that 

From the spherical triangle ODP we have 



whence 



Again 



cos PD — cos $ cos Sv 

cos y = r-^-T-^r— -^ 

sm 6 sm o^ 



sin 6 80 = cos y sin ^ Sx + 2 cos d ^x^ •— ^ cos 6 Bd^ 
= cos y sin 6Sx + ^ cos 6 sin^ y S^^. 

sin B^ sin j 



sin Sx sin (6 + S6) 



? 



whence 



86 = "7— T (8y — cot 6 cos y Sv^). 



Substituting these values of SO, S^ in (46) and equating coefficients of Sx^, we 
obtain 

d^w sin^ 7 d^w sin 2^ d^w ^ d^w mi2^ q,ob6 dw . ^ ^dw 

d^ "" sin^^ #2 + sin elied4'^ ^^^ ^ ^7^ ~^^f6^ I4 + ®^^ ^^^^^^dB ^^^^* 

Whence it follows, that if p^, p^ are the principal radii of curvature along and 
perpendicular to a meridian 



1 
Pi 


1 
a 


1 




1 

Pa 


1 


1 


( 1 (^^W 



+ cot d ~'\' w\ 



<> * • 



• (48). 



^ 
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23. When the middle surface is inextensible, it has been shown by Lord Rayleigh^ 
that the displacements are given by the equations 



to = — tAs e''^ sin tan* l 9 
v= StA, e^'^^sin^tan^l^? 
w = tAs e''"^ (s + cos 6) tan" ^ 6 



> 



(49). 



where 5 = 2, 3, 4 . . . and A^ is a complex function of the time. From these equations 
it can easily be shown by mea,ns of (48) that 

11 .A, (s^ — s) e^"^ tan* ^0 / 1 1\ ^ , 

~ _ 2 —_--- _ :::z: — — — - . . . . (50). 

The value of the potential energy is given by the second line of (16) ; also by the 
first two of (8) and by (48) 

X=^/x=---- (51), 

and by the last of (8) 

___ ^^ ^ A, {s^ — s) 6''^ tan" ^ $ 
^ "^ a sin^ 6 

whence 

^ = Zj^e ^^^ ^' ^'' ~ ') "*"' '"^ ^^""^ ier+{tA, {s^ - s) sin 4 tan^ ^ 9}'] (52), 

which agrees with Lord Ravleigh's result. 

Let us now suppose that a bell which consists of a spherical shell, bounded by a 
small circle whose latitude is ^tt — a, is vibrating in such a manner that its middle 
surface does not undergo any extension or contraction throughout the motion. One 
of the boundary conditions requires that the flexural couple G^ should vanish along 
the circle of latitude which constitutes the free edge of the bell. By (7) and (51) 

G, = inh^m = ^n¥ {\ + E(X + ,x)} 

^ \Pi ^^ 

From (50) we see that G^ cannot vanish for any value of except = 0, that is, 
at the pole, provided s > 2. It, therefore, follows that a spherical bell whose edge is 
free cannot vibrate in this manner if the middle surface is supposed to remain 
absolutely inextensible throughout the motion. If, however, extension or contraction 

^ ' London Mafch. Soc. Proc.,' vol 18, p. 4 (1881). 
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were to take place in the neighbourhood of the edge, it would be possible for Gg to 
vanish there, and also to satisfy the other boundary conditions. 

There seems no reason to doubt that the argument which has been employed in the 
case of a cylindrical shell, would apply equally to the case of a spherical shell, and 
probably also to a shell of any shape ; in which case, the portions of the displace- 
ments upon which extension principally depends, would be small compared with the 
portions upon which bending principally depends, except at points whose distances 
from a free edge are comparable with the thickness. At the same time it would 
be very desirable to obtain the solution of some problem relating to the vibrations 
of a shell whose edges are free, in which no supposition is made as to the relative 
magnitudes of the extension al and flexural terms. 



